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1 Introduction 

There has recently been much interest in the odderon, both from experimental and 
theoretical points of view. On the experimental side, various processes mediated by a 
C-odd object with an intercept around unity are planned or already under investigation 
at HERA: the pseudoscalar meson production in ep collisions offers a direct probe 
for the odderon [1]. On the theoretical side, estimates of the relevant cross-sections 
have been made, using simplest model for the odderon, just a C-odd state of three 
non- interacting gluons [2] or a non-perturbative Regge pole at j — 1 [3]. Parallel 
to this there has been much activity in studying the gluon interaction effects, which 
presumably change the odderon intercept from exactly unity. In the perturbative 
QCD approach the odderon is a C-odd state of three reggeized gluons, which interact 
pairwise with a well-defined potential. The relevant equation is known since long ago 
[4]. It was shown that it is conformal invariant and splits into a pair of equations 
for the holomorphic and antiholomorphic parts of the wave function in exactly the 
same manner as for the pomeron [5]. An operator, later called q^, was also found, 
which, on the one hand, commutes with the odderon Hamiltonian H and, on the other 
hand, has a much simpler form than H, thus opening a way to simplify the solution of 
the odderon problem considerably [6]. Finally, much effort has been put into directly 
relating the operators q% and H to avoid using explicit wave functions in search of the 
odderon intercept [7-10]. Following this latter line of approach R.Janik and J.Wosiek 
(JW) calculated the odderon energy for arbitrary complex eigenvalue of q% in [8]. In 
a later publication [9] they diagonalized the 53 operator and found both the odderon 
intercept and its wave function. If one relates the odderon intercept to the odderon 
"energy" e as Q 

a o (0) = 1 - (9a s /27r)e, (1) 
the result of J W for the ground state is 

e = 0.16478... (2) 

Thus they confirm our old conclusion that the perturbative odderon intercept lies 
slightly below unity [11]. 

1 Different autors use different definitions of the odderon energy and therefore of its relation to the 
intercept. This has to be kept in mind comparing our numerical values with other references. 
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However a discussion about the validity of their procedure relating q% and H is still 
going on [12]. In particular, for the eigenvalue 93 = the result which follows from 
their procedure is contested [12] (see also Sec. 5). 

In view of this dispute a direct calculation of the odderon intercept, which starts 
from its explicit expression in terms of the wave function, aquires certain importance. 
It resolves in a unique manner any ambiguities involved in formal relations between ill- 
defined operators q 3 and H and associated with boundary conditions to be imposed on 
their respective eigenfunctions. Reporting on these direct calculations in this paper, 
we also present two independent variational results, which can be obtained starting 
from the direct expression of the odderon energy in terms of its wave function. They 
present some interest, since the functions determined in a variational procedure may 
be a good approximation to the exact one, which is quite complicated and ill-suited 
for practical calculations. 

Our direct calculations confirm the value (2) found by JW to be the ground state 
odderon energy. Our variational calculations with two different forms of the trial 
function give 



that is, some 30% higher value for the energy. This change in e corresponds to a 2% 
change in the odderon intercept, for realistic values of a s . From our calculations it also 
follows that the value g 3 = does not lead to any physical odderon state. 

The paper is organized as follows. In Sections 2 and 3 we review the basic formulas 
which serve as a starting point for our calculations. Sections 4 and 5 present the results 
of direct calulations of the odderon energy with a variational and exact wave functions 
respectively. Some conclusions are drawn in Section 6. 

2 The odderon energy 

The odderon energy can be sought as a ratio 



where E and D are energy and normalization functionals, quadratic in the odderon 
wave function Z(r,(f>) [13]. Explicitly 



e = 0.223 and e = 0.226 



(3) 



e = E/D 



(4) 




(5) 
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Here 

e n {v) = 2 Re + w) - ^ , (6) 

a n (z/) is a double Fourier transform of Z(r, 0): 

/■oo 

a n (z/) = / drr^ 2w / rf0 e - m * 
Jo Jo 

(■ n+1 jrh/1 . n — l.\,. n — l w ~ . n — l NrW 

+ — + re»*(/i -iv- — )J (w/ - + ^ - —)^(^ ^) (7) 

where /i and /i are the two conformal weights for the holomorphic and antiholomorpic 
parts of the wave function. They have a general form 

h = 1/2 + n — iv, ft, = 1/2 — n + iv, — oo < v < oo, n — . . . — 1, 0, 1, . . . (8) 

For the lowest branch of the odderon spectrum, which leads to the highest intercept, 
supposedly n — 0. This will always be assumed in the following. For the ground state 
also v — 0. D is obtained substituting e n (u) by unity in (5). 

The odderon wave function Z can be considered as a function of z — rexp(i0) and 
its conjugate. It has a form 

Z(z,z*) = \z(l-z)\ 2h/3 $(z,z*) (9) 

where function $ has to be invariant under substitutions 

z^l-z, z^l/z, (10) 

which follows from the requirement of Bose symmetry in the three gluons. 

Eqs. (4), (5) present a direct way of calculating the odderon energy once its wave 
function is known. It also opens up a possibility of finding e by variational techniques 
minimizing the value of (4) in a given space of functions which should satisfy (9) 
and (10). Note that if one drops the restricting condition of Bose symmetry (10), 
then function Z and consequently a n (v) become abritrary, with the only requirement 
that (5) exist. Then it becomes clear that the minimal value of E/D is realized by 
l^n^)! 2 = 8 n oS(v), which gives e = eo(0) = —4 In 2, that is the pomeron energy. Since 
the space of functions obeying (10) is smaller, one gets an evident bound [14] 

e>-41n2 (11) 



In fact this bound is very crude (cf. (2)). 
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3 The (J3 operator 

An operator q 3 which commutes with the odderon Hamiltonian was found in [6]: 

Ql = -rnrlsrlqjqlqi [g 3 2 , H] = 0. (12) 

Evidently, the odderon ground state (nondegenerate) should also be an eigenstate for 
q 3 . But, in contrast to the odderon Hamiltonian H, the operator gf is a finite order 
differential operator, which does not contain logarithms of neither momenta nor co- 
ordinates. It splits into a product of two differential operators q 3 and its conjugate, 
each of the third order in complex variables z and z* respectively. In this section we 
closely follow the derivation of the spectrum of q 3 given in [9] . The eigenvalue equation 
for 93 and its conjugate can easily be obtained using the explicit form (12). For the 
holomorphic part of the function $ (see Eq. (9)) it reads [6,9] 

a(z)^(z) + b(z)j^*(z) + c{z)±*{z) + d(z)*(z) = 0, (13) 

where 

a{z) = z 3 (l - zf, b(z) = 2z\l - zf(l - 2z), 
c(z) = z(z - 1) (z(z - l)(3/i + 2)(/i - 1) + 3/i 2 - /i) , 
d(z) = /2 2 (1 -n){z+ l)(z - 2){2z - 1) - iq 3 z(l - z) 

and \x = h/3. 

This is a third order linear differential equation with three regular singular points 
at z — 0, 1 and 00. It has three linearly independent solutions uf \ i = 1,2,3 which 
can be chosen so as to possess a given behaviour in the vicinity of z = 0: 

uf\z)^z 2h / 3 , u { ?\z)~z l - h '\ u { 3 0) ~ z i-W\ogz + az- h / 3 , z^O (14) 

obtained from the characteristic equation corresponding to (13) at small z. 

The final eigenfunction z*) should be constructed as a sum of products of these 
solutions with those for the antiholomorhic part. 

3 



H>( z , z *) = £ uf\z*)A^uf\z) = u^A^u^ (15) 
;,fc=i 



Apart from the Bose symmetry requirements (10) it has also to be a single-valued 
function in the complex z, z* plane. This latter condition puts evident restrictions 
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on the form of the coupling matrix . From the behaviour (14) one finds that the 
function $ will be uniquely defined in the vicinity of z = if and only if 

4(0) _ A (P) _ AO) _ AO) _ n AO) _ AO) (m] 

A 12 — A 2l — A 13 — ^31 — u ; ^23 — ^32 l i0 J 

To have uniqueness around the two other singular points, 1 and oo, one has to know 
the behaviour of (15) in their vicinity. Due to the symmetry of the Eq. (13) under the 
substitutions 

z->l/(l-z) (17) 

one can easily construct two other systems of solutions u^p and u\°°\ i = 1,2,3 with 
the same behaviour around z — 1 and z = oo respectively. They can be expressed as 
superpositions of the initial solutions,say, 

4 1) W = E^ 0) 4 0) W (is) 
k=i 

and similarly for u^°°\z). The "transfer matrix" from the solutions i^ ) to solu- 

tions is a constant matrix, which can technically be calculated once the solutions 
and are known. Evidently, in terms of the wave function $ can be written 

as 

= u^A {l) u {l) (19) 

where 

A {i) = (ij(oi))T A (o) jR (oi) ) ^(oi) = (i2 (io))-i (20) 

For $ to be a single valued function of in the vicinity of z — 1 it is necessary that 
has the same properties (16) as the matrix A^°\ Moreover the Bose symmetry 
requires that these matrices coincide. This gives an equation 

A = (R<M) T AR<M (21) 

for a matrix A satisfying (16) [9]. It determines both the eigenvalues q$ and non-zero 
elements of the matrix A, that is, the eigenfunction <&(z,z*) according to (15). Note 
that there is no need to additionally require that $ should be a single valued function 
of around z = oo, since a contour encircling this point can be made of two contours 
around z = and z — 1. 

In [9] the transfer matrix R was determined numerically from the solutions u cal- 
culated as a series in powers of z. Solution of (21) then determined the eigenvalues of 
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qs- They are all pure imaginary. The one corresponding to the odderon ground state 
turned out to be 

q 3 = -0.20526 i (22) 

After that the value (2) for the odderon energy was obtained by JW using their pro- 
cedure to relate the eigenvalues of and H constructed in [8]. 

Our calculations use the form of the wave function (15) with the matrix A deter- 
mined by Eq. (21) as an input to be substituted into (5), which gives the value of the 
odderon energy directly. 



4 Variational calculations 

4.1 Starting point [13,14] 

The first calculations of the odderon energy starting from Eqs. (4), (5) were done in 
the variational technique. In fact, taking some trial function which satisfies (10) and 
putting it into (4) and (5) one obtains an upper bound for the odderon energy and 
thus a lower bound for its intercept. Of course the problem consists in choosing a 
good trial function, on the one hand, close enough to the exact one and, on the other 
hand, sufficiently simple to allow for the numerical treatment. On top of that there is 
a problem of satisfying conditions (10). 

This latter problem can be resolved by choosing the trial function as a function of 
arguments which are invariant under (10) by themselves. One of such arguments was 
proposed in [13]. Let r = \z\ and r\ — |1 — z\. Then it is easy to see that 



^2^,2 



(23) 



(l + r 2 )(l +r 2 )(r 2 + r 2 ) 
is invariant under (10). This is not the only argument with this property. Let us take 

, = (i-M 2 )(i-N 2 )(M 2 -N 2 ) m] 

(l + r 2 )(l + r 2 )(r 2 + r 2 ) 1 1 

It is invariant under z — > 1 — z and changes sign under z — > 1/z, so that b 2 is invariant. 
Thus any function $(a, b 2 ) will satisfy (10). To further narrow a possible choice of the 
trial functions, one can impose condition that at z — * it has the correct behaviour 
following from (14) and (15). Due to invariance under (10) it would mean that it would 
also have the correct behaviour near the other two singular points z — 1 and z = oo. 



8 



For the ground state h — 1/2 it means that at z — > 

$(z,z*) ~ C!r 2/3 + c 2 r 5/3 ( 1 + c 3 In r + c 4 cos 20) (25) 
This suggests taking the trial function in the form 

$(a, b 2 ) = Cl a 1/3 + c 2 a 5/6 (l + c 3 In a + c 4 6 2 / a ) (26) 
In our calculations the trial functions were chosen in a more general form 

N—Ni Ni 

$ = £ Cka"' 2 - 1 ' 6 + £ 4a fc " 1/6 In a + /6 2 /a" 1 / 6 (27) 
fc=i fc=i 

with A/" + 1 variational parameters c&, dk and / (one of them is determined by the nor- 
malization condition). Note that the last term in $ contains an azimuthal dependence, 
which, due to the properties of e n (v) , can only raise the energy. Therefore one cannot 
expect any improvement of the variational bound coming from it. Indeed, our calcu- 
lations incuding this term always lead to / close to zero for the minimizing function. 
Accordingly in the following we shall not discuss its influence and use (27) with / = 0. 
Function $ with AT = 3, Aq = 1 was used in the first variational calculations in [13]. 
The basic quantity a given by Eq. (7) can be presented in the form 

«»M = {\ ~ ™)C^ " ^)[(^ + ^V) + ^/SiMl (28) 

where 

fOO c-Itx 

ff)(v)= drr~ l - k - 2lv d(t)e- in4, Z{r )( t)) (29) 
Jo Jo 

Function f!^\v) has the properties: fn\v) = f^(—v) = (f^)*- Using them one 
can restrict the summation over n and integration over v to nonnegative values. The 
value of the |a n (z/)| 2 can evidently be expressed via a single function f^(u), which will 
simply be denoted as f n (v) in the following. In this manner one obtains 

E = ]T / dve n {y)Vn(y) (30) 

where for n > 

VnV) = 4 + ^^T- + V 2 )[(^- + V 2 )\U»)\ 2 + + - 2 )IA-l(-)| 2 + 

71 -\- 1 2 — 71 

2 Re (-^- - iu) (— " ^)/„H/„-lH] (31) 
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and 

Po(^) = ^ 2 (i + ^ 2 )[(i + ^ 2 )|/ (^)| 2 + (l + ^ 2 )|/i(^)| 2 + 2Re(i-^)(l-^)/ (^)/ 1 (^)] 

(32) 

The normalization functional D has the same form (30) with e n (v) — > 1. Thus calcu- 
lation of the odderon energy requires calculation of functions f n {v) and e n {v). 

The main technical difficulty is the double Fourier transform (29). The energy e n (u) 
in E, Eq. (6), monotonously grows both with n and v. It is negative only for n = 
and small enough values of v. So the problem with this formalism is that cutting in 
(5) summation over n and integration over v by some maximal values n m and u m , 
one always gets E smaller than the exact value, corresponding to n m and v m — > oo. 
Therefore in the course of the calculation one always approaches the variational value 
of e from below. As we shall see in the following, in fact, rather high values of n m and 
v m are necessary to obtain e with a good degree of accuracy. This is the reason why 
in [13] a negative value was obtained for e (corresponding to «o(0) ~ 1.07 ): too small 
values of n m and u m were chosen there. On the other hand, with high n and v, the 
double Fourier transform (29) becomes very difficult, especially having in mind that, 
due to the factors in (31), two first terms in the asymptotic expansion of a at high n 
and v cancel. As a result, a trustworthy calculation of E and D turns out to be very 
complicated, in spite of its superficial transparency. 

The important point in obtaining reliable results has been using analytic asymptotic 
expansions for / at high n and u, which are briefly discussed in Appendix. 

4.2 Numerical procedure 

As mentioned numerical integration in the double Fourier transform (29) presents a 
formidable calculational task. As mentioned, our results were obtained by two groups 
working independently and using different choices of the number JVi of logarithmic 
terms in (27) for a given total number of terms N. We present here in some detail the 
calculational procedure adopted in the computation with Ni — 1. 

In the integral (29) the integration over r > 1 was transformed to r < 1 by a 
substitution r — > 1/r. The integration over <fi was reduced to the interval < <fi < n 
and exp(— n<f>) was substituted by 2cos(n0). To soften the behaviour of the integrand 
at small values of r, three first terms of its asymptotics at r — > were subtracted and 
treated in an exact manner. The final formula for f n (v) used in the calculations was 
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thus 

/ n (i/) = 2 f #cosm0 /'^(^(C^^-Co) 

JO J 

+r- 1+2 -(C(r, 0) - Co + Cir In r - C 2 r)) + 

^ 2 S«/(i-2w/) + (i + r+k; (33) 

where ( = Z/r and the subtraction constants are 

Co = c 1 (l/2) 1 / 3 , Ci = -2c w (l/2) 5 / 6 , C 2 = (l/2) 5/6 (c 2 - Cjv ln2) (34) 

Eq. (33) was used for numerical calculation of f n {v) in the interval of < n < 30 and 
< v < 15. Integrations were performed by dividing the rectangle < r < 1, < 
< 7r into an M x M grid, interpolating C quadratically on the grid and then doing 
the integrals explicitly. The maximal value of M was 640. The achieved accuracy 
was about 10~ 5 . Thus calculated values of f n (v) were summed over n and integrated 
over v as indicated in (30)-(32) to obtain E and D. Stable results were obtained with 
the quite high maximal values n m = 300 and u m = 150. In the part of (n, v) space 
outside the rectangle < n < 30, < v < 15 the asymptotic expressions were used 
for f n {v) (see Appendix). As a result we calculated E and D as a quadratic form 
in the variational parameters Ck,dk- Afterwards the minimal value e of E, subject to 
condition D — 1, was found by standard methods. 

Other calculations used (27) with = 1 for TV = 3, 4 and N ± = 2 for N = 6. 
The adopted numerical procedure was different, however the final results, as we shall 
presently see, are quite similar. 

Our results for different number of parameters iV are presented in Table 1 for the 
both choices: Aq = 1 always and Aq = 1 for N = 3,4, Aq = 2 for N = 6. The 
corresponding energies are denoted e\ and e 2 respectively. Adjacent columns present 
values of the variational parameters for both cases (ci = 1). The standard precision 
corresponds to the (r, 0) grid 320x320. To clarify the accuracy achieved we also present 
the results for e\ with a double precision (the grid 640x640) for A^ = 5. 

Inspecting these results we see that the final accuracy in energy is of the order 
5.10~ 3 . Also taking A^ > 6 with Aq = 1 leads to no improvement within the precision 
achieved, since the corresponding change in energy is of the same order or less. So our 
conclusion is that the variational odderon energy with a trial function (27) is given 
by (3) and that with the accuracy achieved in the course of numerical integration, as 
described above, the maximal number of terms to be taken in the trial function is 
N = 6, although already with A^ = 3 the energy is obtained up to 1%. 
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5 Exact wave functions 
5.1 Problems of precision 

In this section we report on the calculations of the odderon energy using Eqs. (4) 
and (5) with the odderon wave function (15) obtained after solving Eq. (21). These 
functions can be obtained in the form of power series in z and z*, different in different 
parts of the complex z, z* plane. The exact solution is, of course, analytic in each of the 
variables, except at the three mentioned singular points z — 0, 1, oo. So it is absolutely 
smooth at the boundaries of the regions in which it is represented by different power 
series. However, in practice one knows it to some finite precision. This causes certain 
discontinuities in the function itself and its derivatives at the mentioned boundaries. 
At high values of n and v the double Fourier transform (29) is very sensitive to such 
local irregularities of the wave function. They result in a very poor precision for the 
calculated odderon energy and even in a completely wrong order of magnitude for it. 
Because of this, our first step was to obtain the odderon wave function with a higher 
precision than reported in [9]. To this end we set up a program which essentially repeats 
the procedure employed in [9] and allows to reduce discontinuities of $ calculated in 
different variables to values of the order 10~ 9 . The only difference is that we used 
the basic functions Ui, i — 1,2, 3, multiplied by appropriate factors to make them 
real at points where the transfer matrices are calculated. This substantially facilitates 
achieving the desired accuracy. 

For the eigenvalue and wave function parameters a, f3 and 7 we obtained the 
following values (precision 10~ 9 ) 

iq 3 = 0.205257506, a = 0.709605410, (3 = -0.689380668, 7 = 0.145651837 (35) 

However even with these high precision values direct calculation of $ in one of the three 
sets of variables z, 1 — l/zorl/(l — z) fails in the vicinity of the point zo = exp(i7r/3) 
where none of the series converges absolutely. In spite of the fact that zq is only a 
single point in the (r, 0) plane, this makes it practically impossible to calculate the 
double Fourier transform for |n| > 10 and/or \v\ > 5. To overcome this difficulty we 
had to redevelop the function $ around the intermediate point (1/2)zq. With this 
redevelopment 100 terms in the series proved to be sufficient to obtain reliable results. 

The Fourier transform was again performed by interpolating Z quadratically on 
the (r, 0) grid and doing the integrals in r and analytically. Reasonable results are 
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obtained already with a 160 x 160 grid. We however also used 320 x 320 and 640 x 640 
grids to analyse the precision achieved at this step. 

Even with a 640 x 640 grid the numerical Fourier transform becomes unreliable 
for | n | > 30 and/or \u\ > 15. For such high values of |n| and \v\ we used asymptotic 
formulas for the Fourier transform, which can easily be obtained from the expansion of 
$ around z = 1 in a similar manner as for variational wave functions in the preceding 
section (see Appendix). 

Our final cutoffs were chosen to be |n| < 300 and \v\ < 150, which proved to be 
quite sufficient for the determination of e with a precision 0.001. 

5.2 Ground state energy 

Results of our numerical calculation of D and E in the region \n\ < 30 and \v\ < 15 
using an N x N grid in the (r, 0) plane are presented in the Table 2 for N = 160, 320 
and 640. 

To these values one has to add the contributions from the asymptotic region de- 
scribed in the preceding section. They are 

AD = 0.002398, AE = 0.018451 (36) 
Taking the results at iV = 640 as the most accurate ones we finally have 

D = 1.644454, E = 0.273083, e = 0.1660 (37) 

Thus our result for e coincides with the value found in [9] up to 10" 3 . With all the 
difficulties involved in the numerical calculations, we consider this agreement quite 
satisfactory. So direct calculation of the odderon energy confirms the result found by 
JW. 

5.3 Higher eigenvalues of q% 

We have tried to check the result of [9] also for the excited state with the next higher 
value of iq%. Unfortunately in this case calculations proved to be still more difficult 
and we could not arrive at a result of a convincing accuracy. 

Our precise calculations of the wave function gave for this state: 

iq 3 = 2.343921063, a = 0.391855163, (3 = -0.0533712012, 7 = 0.918477570 (38) 
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Numerical calculation of D and E in the region \n\ < 15, \v\ < 15 gave results presented 
in Table 3. As one observes, the achieved accuracy does not exceed 15%. Analyzing 
these numbers one can see that all the error comes from the region of maximal \n\ and 
\u\ where the double Fourier transform is apparently performed inaccurately. ^From 
these numbers we can only conclude that for this excited state 

e~ 2.0 ±0.3 (39) 

The value found in [9] is 1.7231... Our result does not contradict this number. 

5.4 The case q s = 

We have also studied a degenerate solution for q 3 = 0, discussed in [9]. However in 
this case we were not able to construct a wave function unique in the (r, 0) plane and 
satisfying the necessary boundary conditions and symmetry requirements. At q 3 = 
the solution uf\z) has no logarithmic term. If, following [9], we seek $ in the form 
(15) and require the same conditions (16) to be satisfied, then Eq. (21) has no solutions 
at all, which means that the Bose symmetry cannot be fulfilled. The problem is related 
to the fact, that, with the logarithmic term missing in u 3 , a unique solution is obtained 
also with A 33 7^ 0. As a result, starting from a solution with A 33 = one always 
obtains A 33 7^ after applying the transformation (21). 

One may wonder if a solution containing a product u 3 (z*)u 3 (z) is admissible as 
a physical odderon wave function. Such a solution does not vanish as the distance 
between any of the three gluons becomes small. This can be easily seen from the 
behaviour at z — > of the complete wave function (9): the explicit factor r 2/l//3 cancels 
against the leading term r ~ 2h ^ 3 in the product u 3 (z*)u 3 (z) and one gets a constant. In 
fact the solution with this property can be easily constructed explicitly [12]. It is just 
a sum of three pomeron wave functions 

r 2 , r 3 ) = ^ B FKL(r 1 ,r 2 ) + ^bfkl{t 2 , r 3 ) + ^BFKi(r 3 , n) (40) 

This function is an evident eigenfunction of q 3 for zero eigenvalue and it is conformal 
and Bose symmetric. However it cannot be considered as a physical odderon state, 
since the Hamiltonian H cannot be applied to it due to its singularities at low gluon 
momenta. Take the first term, which does not depend on r 3 . Then it is proportional to 
5 2 {k 3 ) in the momentum space hi, k 2 , k 3 of the three gluons. The part of the odderon 
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Hamiltonian depending on k% contains a term proportional to In k\ and two interactions 
between the gluons 31 and 32. It can be easily seen that the latter give finite result, 
applied to a wave function containing 5 2 {k^). The ln&f, however, is infinite. So H 
applied to (40) is not defined. This property can also be seen from the approximate 
relation H ~ lng 3 valid for small values of z [5]. Evidently H diverges as g 3 — > 0. 

Thus our conclusion is that the qs = state does not correspond to any physical 
odderon state. 

5.5 "Moving" odderon 

For conformal weights h = \ + ia the odderon energy is supposed to behave at small 
a as 

e(cr) = e + aa 2 (41) 

where e is the value (2) and a is a parameter which determines the diffusion of the 
odderon wave function in the momentum space. This parameter has long been known 
for the pomeron to be 14£(3) ~ 16.8 (in units 3a s /n). It is of certain interest to find 
a for the odderon. To this aim we first found the parameters of the odderon wave 
function for various (small) o using the same method as employed for h = 1/2. Our 
results are presented in Table 4. The value of iq^ turned out to be real for arbitrary a, 
whereas, with a chosen to be real, both (5 and 7 result complex. We chose a — 1. 

Inspecting these figures one immediately notes that I7I = iq 3 . This relation was 
predicted by L.N.Lipatov [15]. 

With the odderon parameters found, we calculated the odderon energies directly, 
using the same techique as for a = 0. With a different from zero calculation becomes 
still more cumbersome and time and memory consuming due to lack of certain sym- 
metries and overall complex arithmetics. For these reasons we had to limit ourselves 
with a maximal 160 x 160 grid in the (r, 0) plane and neglected the contribution from 
the asymptotical region n > 30, \u\ > 15. Our results are shown in Table 5 together 
with the ones obtained via the solution of the Baxter equation [16]. 

Our energies lie a little below the ones obtained from the Baxter equation, which is 
natural since we have neglected the asymptotic part of the n, v. Having this in mind 
we find a complete agreement between our direct calculation results and the ones based 
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on the Baxter equation. ^From our energies we find for the parameter a in (41) 

a = 0.61. 

More precise energies found in [14] lead to 

a = 0.605. 

Note however that already at a — 1 the approximation (41) breaks down and more 
powers of a 2 are needed to describe the energy behaviour. It is interesting that the 
parameter a for the odderon is much smaller than the one for the pomeron. In fact 
their ratio is of the same order as the ratio of corresponding energies. 

6 Conclusions 

Our calculations confirm the results obtained by JW for the perturbative odderon 
intercept and thus seem to remove any doubts concerning the validity of their procedure 
to relate the operator and the Hamiltonian. 

The variational calculations give a result for e which is ~30% larger than the exact 
value. However they also convey the important message that the intercept of the 
odderon lies quite close to unity being slightly smaller. In fact, for realistic values of 
a s , the intercept ao(0) calculated variationally is at most 2% lower than the exact one: 
0.94 instead of 0.96. 

The disputed eigenvalue g 3 = does not seem to correspond to any physical odderon 
state. 

In the course of the variational calculations a simple approximate form of the odd- 
eron wave function is found, which allows a realistic calculation of the odderon residues, 
important for the study of processes involving the odderon exchange. 
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8 Appendix. Asymptotics at large n and v for vari- 
ational calculations 

Passing to the variable p = — In r and introducing 2-dimensional vectors x = (p, 0) 
and w = (z, n) = (2v>, n) we rewrite (29) as 

dp / d<f>e iwx Z(x) (42) 

-OO J — IT 

The integration point x = is obviously essential for the asymptotics at high n and v. 
At x — > 0, keeping terms up to third order in small p and 0, we have (x = \J p 1 + 2 ): 

Z p = (rrx) 1 / 3 ^ = (l/2)V p+1 / 3 (l - (l/2)p + a lP 2 + b^ 2 + c lP 3 + d lP <p 2 ) (43) 

where 

ai = 5/36 - (29/12)p , 6i = -1/72 - (25/12)p , 
Cl = -1/36 + (29/24)p , di = 1/144 + (25/24)p 

and 

p = — 1/6 , A; = 1,2,... 
For the term with a logarithm, in the same manner we obtain 

Z p = (rrx^V In a = (l/2)V p+1 / 3 [ln(x 2 )(l - (l/2)p + a lP 2 + b^ 2 + c lP 3 + rf x p0 2 ) + 
ln(l/2) - (1/2) ln(l/2) p + a 2 p 2 + 6 2 2 + c 2 p 3 + d 2 p<p 2 )} (44) 

where 

p= fc- 1/6 , k = 1,2,... 

and 

a 2 = ^ ln(l/2) -29/12 , fo 2 = 6i ln(l/2) - 25/12 , 

c 2 = ciln(l/2) + 29/24 , d 2 = di ln(l/2) + 25/24 

Inserting these expressions into the integral (42) and extending the integration over 
to the whole real axis one obtains the asymptotical expansion of different terms in 
fniy). In particular the asymptotical expansion of the term originating from Z p is 
found as 

f<p\v) = ^(i + (l/2)id/dz - ad 2 /dz 2 - M 2 /dn 2 + icd 3 /dz 3 + idd 3 / 'dzdn^w^- 1 
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where a = 2p + 4/3 and 

Cp = 2 a+1 - p T 2 (l/2 + a/2) cos(ira/2) 

Doing the derivatives, one obtains finally (w = \Jn 2 + Au 2 ): 

Re/^(z/) = CpW""- 1 (l + (a + l)w"V( & - a(a + 2) + n 2 (a - b(a + 2))]) (45) 

Im/^(z/) = -(l/2)c p (« + 1)^ Q - 3 
(l + 2(a + 3)ur 4 [z 2 ( c (" + 2) - d) + n 2 (d(a + 4) - 3c)]) (46) 

For the term with a logarithm only the part with lnx 2 contributes. The result 
coincides with the formula above where the constant c p is substituted by 167r(l/2) p .If 
one puts these asymptotic expressions into (31) one finds that the two leading terms 
coming from Zi/ 3 cancel. Numerically the asymptotic expansion begins to work at 
rather high values of n and v: \/n 2 + 4z/ 2 > 30. 
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Table 1. Odderon energy and parameters of the trial 

functions 



N 


ei (N x = 1) 


c 2 , c 3 , di 


e 2 (iVi = l,2) 


c 2 , c 3 , ...; di, d 2 


3 


0.22865 


-0.5036; 0.2895 


0.23137 


-0.50276;0.28936 


4 


0.22632 


-0.2791,-0.3190;0.3609 


0.22735 


-0.21420,-0.41031;0.38136 


5 


0.22627 


-0.2009,-0.5052,0.1557; 
0.3779 






5* 


0.22634 


-0.2021,-0.5028,0.1543; 
0.3775 






6 


0.22619 


-0.3735,0.08842,-0.9765, 
1.003;0.3471 


0.22269 


0.49231,-3.49272,-1.82821; 
0.50316,-1.49528 


7 


0.22618 








8 


0.22616 








9 


0.22616 









*) double precision. 



Table 2. D and E for the ground state 



N 


D 


E 


160 


1.642162 


0.255480 


320 


1.642085 


0.254852 


640 


1.642056 


0.254632 



Table 3. D and E for the state with iq 3 = 2.343921063 



N 


D 


E 


160 


2.92863 


6.08989 


320 


2.80693 


5.27381 


640 


3.05939 


6.96764 



Table 4. Odderon parameters for h — \ + ia 



a 






7 


0.01 


0.205306079 


-0.971740164-i0.014404102 


0.205305637-i0.000425478 


0.1 


0.210089247 


-0.995153863-i0. 142974530 


0.210052319-i0.003938872 


0.3 


0.247227544 


-1.156524786-10.415163678 


0.247186043-i0.004529717 


0.5 


0.316528176 


-1.395571390-10.695891904 


0.316214188+i0.014095104 


1.0 


0.619239545 


-2.044631201-il.672240784 


0.591391973+i0.183611401 



Table 5. Odderon energies for h — \ + ia 



a 


e 


6 [14] 


0.0 


0.1534 


0.16478 


0.1 


0.1597 




0.3 


0.2085 


0.21777 


0.5 


0.2980 


0.30523 


1.0 


0.6269 


0.63228 



